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Abstract 

In this paper, we give a direct method for calculating the partition function, and hence the 
equation of state (EOS) of Quantum Chromodynamics (QCD) at finite chemical potential and 
zero temperature. In the EOS derived in this paper the pressure density is the sum of two terms: 
the first term V{^)\^=q (the pressure density at /x = 0) is a //-independent constant; the second 
term, which is totally determined by G r [fj] {p) (the renormalized dressed quark propagator at finite 
/i), contains all the nontrivial //-dependence. By applying a general result in the rainbow-ladder 
approximation of the Dyson- Schwinger approach obtained in our previous study [Phys. Rev. C 71, 
015205 (2005)], G'i?[//](p) is calculated from the meromorphic quark propagator proposed in [Phys. 
Rev. D 70, 014014 (2004)]. From this the full analytic expression of the EOS of QCD at finite // 
and zero T is obtained (apart from the constant term 'P(/i)|^=o which can in principle be calculated 
from the CJT effective action). A comparison between our EOS and the cold, perturbative EOS 
of QCD of Fraga, Pisarski and Schaffner-Bielich is made. It is expected that our EOS can provide 
a possible new approach for the study of neutron stars. 
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The study of the partition function is at the crux of equihbrium statistical field theory 
(see, for example, Refs. [1,2]). The thermodynamic properties of a system, hence also 
the equation of state (EOS), are completely determined by the partition function. The 
calculation of the partition function of Quantum Chromodynamics (QCD) at finite chemical 
potential is a contemporary focus; e.g., Refs. [3-9]. In addition, it is well-known that in 
astrophysics the study of the neutron star depends crucially on the assumed EOS [10,11]. 
The study of EOS of QCD is thus of extreme importance. In this paper, we try to give 
a direct method for calculating the partition function and EOS of QCD at finite chemical 
potential. 

The renormalized partition function of QCD at zero temperature and finite chemical 
potential reads 

Z[^x] = j Vq-RVgnVAR expl^-SR[qR,qR,AR] + J d^x //Z2^i?(x)74gj?(x)| , (1) 

where SR[qR,qR,AR] is the standard renormalized Euclidean QCD action with qR being 
the renormalized quark field with three flavors and three colors, Z2 — Z2(C^,A^) is the 
quark wave-function renormalization constant (C is the renormalization point and A is the 
regularization mass-scale). Here we leave the ghost field term and its integration measure 
to be understood. The pressure density V{fi) is given by 

7'(/i) = i lnZ[/.], (2) 

where V is the four-volume normalising factor. The above equation for the pressure is just 
the EOS and from this one immediately obtains the quark-number density 
dV{n) 1 1 dZ[n] 



^^^^ dfi VZ[fj] dii 

1 J VqRVqRVAR J d^xZ2qR{x)-f4qR{x) exp {-SnlqR, qR, Ar; ij]} 



V JVqRVqRVAR exp {-SB\qB,, qn, Ar- fi]} 

where SrI^r, qR, Ar; fj] = SrI^r, qR, Ar] - J d'^x lJ,Z2qR{x)'y4qR{x) . 

On the other hand, the dressed quark propagator at finite chemical potential can be 
written as 

^ r 1/ X _ JVqRVqRDAR qRi{x)qRj{y) exp {-SR[qR,qR, Ar; /j,]} 
UR,,inl[x,y) - JVqRVqRVAR exp{-SR[qR,qR,AR;n])} " 

Prom Eq. (4) , it is easy to obtain the following 

rT.fri!^ 1 I'^Qr'I^Qr'^ArJ d'^xqR{x)-f4qR{x)cxp{-SR[qR,qR,AR;^]} 

Ir {<^Rmi4 = r-n- t. /i r c r- ~a n ' 

j VqRVqRVAR exp {-SR[qR, qR, Ar; /ijj 
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where the notation Tr denotes trace over the color, flavor, Dirac and coordinate space indices. 
Comparing Eq. (3) with (5), we obtain a well-known result (for its recent application, see, 
e.g. Ref. [12]) 

= -Y^r{Gn[M = -N^NfZ, J ^tr {G^[/.](p)74} , (6) 

where Nc and Nf denote the number of colors and of flavors, respectively, and the trace 
operation is over Dirac indices. Prom Eq. (6) it can be seen that the quark-number density 
p{li) is totally determined by the dressed quark propagator at flnite chemical potential. 
Setting /X = in Eq. (6), one flnds that the quark- number density at zero chemical potential 
vanishes. This is because by writing out the general Lorentz structure of GR[fi = 0]{p) and 
performing the trace one can verify that the integrand is an odd function of p4 and therefore 
the integration vanishes. This is what one expects in advance. 
Integrating the equation = ^^^^ oiie obtains 

-Pifi) = P(//)|^=o + = r{fi)\,=o-NcNfZ, j'di,' j ^tr{G«[/.'](p)74}. (7) 

From the above equation it can be seen that the pressure density is the sum of two 

terms: the flrst term P(//)|^=o (the pressure density at zero //) is only a //-independent con- 
stant; the second term, which is totally determined by Gr[^]{p), contains all the nontrivial 
//-dependence. Here we note that formula (7) is formally model- independent. However, at 
present it is very difficult to calculate 'P(/i)|^=o and Gr[ii\{p) from flrst principles of QCD. 
So when one uses formula (7) to calculate the EOS of QCD, one has to resort to various 
nonperturbative QCD models. 

Over the past few years, considerable progress has been made in the framework of 
the rainbow-ladder approximation of the Dyson- Schwinger (DS) approach [13-16], which 
provides a successful description of various nonperturbative aspects of strong interaction 
physics. We naturally expect that it might be a useful nonperturbative approach in the 
study of EOS of QCD at flnite chemical potential. In this paper we shall employ this 
approach. 

Now let us turn to the actual calculation of Gr[^]{p). Here we apply the following general 
result proved in Refs. [17,18]: Under the rainbow approximation of the Dyson-Schwinger 
equation (DSE), if one ignores the ji dependence of the dressed gluon propagator (this is a 
commonly used approximation in calculating the dressed quark propagator at flnite chemical 
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potential [14,17-23]) and assumes that the dressed quark propagator at finite /i is analytic 
in the neighborhood oi /i — 0, then the inverse dressed quark propagator at finite chemical 
potential can be obtained from the one at zero chemical potential by the following simple 
substitution [17,18]: 



where p = (p,P4+^A*) and G^^(p) — ij-pA{p^)+B(p^) is the inverse dressed quark propagator 
at II — 0. Here one may ask whether the quark propagator given by the above equation is 
valid for any fi. This question has two levels, those of mathematics and physics. On the 
mathematical level, it should be noted that in deriving Eq. (8) we have assumed that \fi\ 
is smaller than the radius of convergence of Taylor expansion of Gj(^[fj]{p) around fi = 
(for details, see Ref. [17]). However, Eq. (8) holds in the whole domain of analyticity of 
^R^ilAip) complex /x-plane, not only within the circle of convergence of /i expansion. 

This is a result of a well-known theorem in complex analysis [24]: Suppose each of two 
functions f{z) and g{z) is analj^ic in a common domain D. If f{z) and g{z) coincide in 
some subportion D' C D, then f{z) = g{z) everywhere in D. On the physical level, it is 
generally believed that at finite /i there will be a chiral transition. In particular, in the chiral 
limit, this is a first order phase transition, and hence the quark mass function derived from 
the quark propagator will drop discontinuously to zero. This feature is not reproduced by 
the quark propagator in Eq. (8). Hence it should be stressed that this not happening is a 
pure assumption of the model employed in this paper. 

According to Eq. (8), once the dressed quark propagator at yU = is known, one can 
obtain the dressed quark propagator at finite fi by means of Eq. (8). Therefore, in order 
to calculate Gr[ii]{p) using Eq. (8), one needs to specify the form of the dressed quark 
propagator at zero chemical potential in advance. In Ref. [25] , guided by the solution of the 
coupled set of DSEs for the ghost, gluon and quark propagator in the Landau gauge, the 
following meromorphic form of the dressed quark propagator is proposed: 



with rrij — aj + ibj. The propagator of this form has np pairs of complex conjugate poles 
located at aj ± ibj. When some bj is set to zero, the pair of complex conjugate poles 
degenerates to a real pole. The residues rj are real (note that a similar meromorphic form of 



(8) 
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the quark propagator was previously proposed in Ref . [26] , in which the residues in the two 
additive terms are complex conjugate of each other). In the chiral limit, the requirement 
that the dressed quark propagator reduces to the free one in the large momentum limit 
entails that 

np -j^ np 

In this paper, following Ref. [25], we set the renormalization point to be = 16 GeV^. 
Here it should be noted that the quark propagator (9) is obtained from a calculation going 
significantly beyond the rainbow approximation (for details, see Ref. [25]) and not necessar- 
ily a solution to the rainbow DSE for the quark propagator at /i — 0. However, it would be 
valid to assume for the moment that the quark propagator obtained by substituting (9) into 
Eq. (8) is an acceptable approximation for the quark propagator at finite chemical potential 
and see what physical results it will yield. With this in mind, from the form of Gr{p) given 
in Eq. (9) we obtain 



np 
np 



• + rrij i + m*j 



...,^A^)g|^^i^.^$^). (10, 

Substituting the above equation into Eq. (6) and performing the trace, one obtains 



(27r)4 fr[{^J p^ + {P4 + i/J-f + mj ^ J p^ + {p4 + i/i)^ + mf j ' 
Now we need to evaluate the integrals in Eq. (11). We first consider the first integral on 



the right-hand-side of (11) and write 

+00 

' 2 ' 



+C30 

j ^4^ P4 + = [ dp ( dp^ ^ 

J P^ + {P4 + ^A*)^ + m3 J J + (P4 + ^A*)^ + 'm^j 



The integral over p4 can be written as a line integral from — oo -|- i/i to -|-oo -|- i/i in the 
complex plane: 

+00 H-oo+i^ 

^ Pi + ^ f dz- ^ 

+ (P4 + ^A*)^ + ^ z"^ + p^ + m- 

■' — oo+i/i 



2 I ^ _L ^2- 
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The latter integral can be evaluated by means of contour integration. For this purpose we 



choose the following contour (see Fig. 1). The function f{z) 



^ has two poles in 



the complex plane which are located at 



± 



±i 



sgn{-pj)^ 



where m| 



aj + Pji and sgn{—Pj) denotes the sign of —Pj. Among the two poles, one lies 



in the upper half plane, the other lies in the lower half plane. The imaginary part of the 
pole in the upper half plane is 



iiP) = \ 



(f 



rl2 



Oil 



Here one should distinguish two cases: 

(i) /I < u!j{p). In this case, since there are no poles inside the contour (see Fig. 1), the 
integral along the contour vanishes. The integral along the straight segment from —R to it! 
vanishes because f{z) is an odd function. In the — > oo limit, the integral along the two 
vertical segments gives zero contribution because f{z) vanishes at infinity. So in the R ^ oo 

—R+ifj, 

limit, the contour integral has its contribution only from / dzf{z). Therefore one has 



j dzf{z) = 0, fx< Ujip). 



(12) 



-oo+ifj, 



(ii) jJL > ojj{p). In this case, there is one pole inside the contour (see Fig. 2). According to 
the above explanation and applying the residue theorem, one has 

.Ujip) .Zj 



-ii: 



o 



Jiez 



R 



FIG. 1: The integration contour in the complex z plane 
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y dzf{z) = 2m Res(/(2;), z = zj) 



Eqs. (12) and (13) can be combined into the following form 

+ OD+ifl 



j dp4 



-oo ■' — oo+i/i 



j dzf{z) = -7rie{iJ, - ojj{p)). 



(13) 



(14) 



Now, when /i < 



— -, one has 0Uj[p) > /J^jvp, so m this case 



9{t,-u;j{p)) = 0. 



When/.>^/ "-+V^ 



2 , ujj{p) < n if and only if |p| < (//^ ~ i;^ ~ "^j)^ > in this case 



e{ix-ujj{p)) = < 



1, when IpI < (/x^ 



^2 _ ^ _ ^^.^1/2 

0, when |p1 > (// - 



.2 /^l 



a,)V^. 



Therefore one has 



d^p 



Pa + ^/i 



+ (p4 + ijJ^y + "T'j 



0, 



— 7ri dp ^(// — ujj{p)) 
when A* < Y — ^ 



_4^(^2_^_^^.)3/2^ when/.>^/ "^'+v^^ 



47r2i 



4// 



Imz 




Rez 



(15) 



FIG. 2: The integration contour in the complex z plane 
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The second integral on the right-hand-side of Eq. (11) can be obtained by making the 



replacement mj — aj + 
Pa + in 



rrij — aj — (5ji in Eq. (15) and it is equal to the first one: 



+ (P4 + il^Y + rnf 
Substituting Eqs. (15) and (16) into Eq. (11) gives 



3/2 



(16) 



2NcNf 



a 4 



(17) 



V 

In the numerical calculations in this paper, we use three sets of parameters given in Ref. [25], 
which represent three forms of the propagator: three real poles (3R), two pairs of complex 
conjugate poles (2CC) and one real pole and one pair of complex conjugate poles (IRICC). 
These parameters are listed in Table I. 



Table 



I. The parameters used in this paper. They are taken directly from Table II of Ref. [25] . 



Parametrization 


ri 


ai (GeV) 


bi (GeV) 


r2 


as (GeV) 


62 (GeV) 


rg as (GeV) 


2CC 


0.360 


0.351 


0.08 


0.140 


-0.899 


0.463 




IRICC 


0.354 


0.377 




0.146 


-0.91 


0.45 




3R 


0.365 


0.341 




1.2 


-1.31 




-1.06 -1.40 



The dependence of p{fj,) on fj, is plotted in Fig. 3. Just as shown in Fig. 3, the obtained 
quark-number density distribution differs significantly from the Fermi distribution of the free 
quark theory. Physically this is a consequence of dynamical chiral symmetry breaking and 
confinement in the low energy region. We note that when fi is smaller than a critical value 
Ho {hq = 351 MeV, 377 MeV and 341 MeV for the 2CC, IRICC and 3R parametrization, 
respectively), the quark- number density vanishes identically. Namely, // = //q is a singular- 
ity which separates two regions with different quark number densities. This result agrees 
qualitatively with the general conclusion of Ref. [27]. In that reference, based on a universal 
argument, it is pointed out that the existence of some singularity at the point /x = /iq and 
T = is a robust and model-independent prediction. The numerical value of the critical 
chemical potential in pure QCD (i.e., with electromagnetic interactions being switched off) 
is estimated to be {niN — 16 MeV)/A^c = 307 MeV (where niN is the nucleon mass and 
A^c = 3 is the number of colours). The value of the critical chemical potential obtained in 
this study is almost of the same order of magnitude as the estimate in that reference. This 
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FIG. 3: Relation between the quark-number density and the chemical potential 



difference can be attributed to the choice of the parameters of the model quark propagator 
employed in this paper. In fact, as was pointed out in Ref. [25], the choice of the parameters 
of the model quark propagator (9) has some arbitrariness. We expect that our work can 
give further constraints on the model parameters. 



Now let us calculate J d/i'pdi'). From Eq. (17) one obtains 





Jo 



1 / dpi'ipi' 



3^ 

'2 _ « .)3/2 



12 



37r2 







-■1 


2 



'■31 ^d) 



(18) 



where aj, (3j) is: 



3/2 



3(^1 - /^) . 
16 





_ ^^/2 _ 




^l,2_^./2 + ^a] + (3]/2-y 


j^2 _ «./2 _ 





3a A 6 A 
H arctan 



(19) 



Now let us turn to the calculation of P(/i)|^=o- The rainbow-ladder approximation of 
DSE is the stationary point equation for the CJT effective action [28] which, evaluated at 
this stationary point, is [29] 



= -2N,Nf 

where G{p) 



In 



A2(p2)p2 + 52(^2y 



In 



+ l+pV,(p2)}, (20) 



ij ■ pcJvip^) + o"s(p^) is the unrenormalized dressed quark 



i7-pA(p2)+B(p2) 

propagator at = 0. From the model quark propagator (9) one obtains 



np 



+ 



p2 + mf' ' ' 



+ 



rjTTij 

p2 j^*2 



(21) 

Prom Eqs. (20,21) and the parameters of the model quark propagator (9) hsted in Table I, 
one can calculate V{^)\fj,=o. 

The determination of the EOS of QCD is a longstanding problem in strong interaction 
physics. Lattice QCD calculations and phenomenological models try to pin down a useable 
EOS since two decades. It is interesting to compare the EOS obtained in this paper to the 
EOS of QCD proposed in previous studies. Here we shall take one prominent example, the 
cold, perturbative EOS of QCD proposed by Fraga, Pisarski and Schaffner-Bielich in Ref. 



[30]. The pressure density to second order in ag in the MS scheme obtained in Ref. [30] is 
quoted as follows 



Vfps{ij) 



1-2©- G + iV,,„^+(n-HiV,),n^(^)n, (22) 



where G = Gq - 0.536A^/ + Nf In Nf, Go = 10.374 ± 0.13 and A is the renormalization 
subtraction point. The scale dependence of the strong coupling constant as{A.) is taken as 



a,(A) 



2Aln(M) API 
(31 u 



1\2 /32/3o 



5 
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where u = ln(AVA|^), Po ^ 11 - 2Nf/3, /3i = 51 - 19Nf/3, and P2 = 2857 - 5033A^//9 + 
325A^|/27. For Nf = 3, = 365 MeV. The only freedom in the model of Ref. [30] 

is the choice of the ratio A/ fi, which is taken to be 2 in that reference. The perturbative 
EOS (22) is applicable only in the chirally symmetric phase, when the chemical potential 
fi is larger than fi^, the chiral phase transition point. A plot of our EOS and the EOS of 
Fraga, Pisarski and Schaffner-Bielich is given in Fig. 4. Here we note that when applying 
our EOS to the study of neutron star, owing to the boundary condition imposed on the 
surface of neutron star, the constant term 7'(/i)|^=o does not contribute to the mass-radius 
relation, so when comparing our EOS with that of Fraga, Pisarski and Schaffner-Bielich in 
Fig. 4, we do not consider this term. From Fig. 4 it can be seen that in the region of fj, 
studied, the pressure density in the EOS of Fraga, Pisarski and Schaffner-Bielich and our 
EOSs (for the three parametrization of the model quark propagator) increases monotonically 
as fjL increases. In large fj, region the EOS of Fraga, Pisarski and Schaffner-Bielich and our 
EOSs show quahtatively similar behaviors: all of them tend to the free quark gas result 
as fjL tends to infinity. Compared with the EOS of Fraga, Pisarski and Schaffner-Bielich, 
our EOSs tend somewhat more rapidly to the free quark gas result. When /x is less than 
about 1.8 GeV, the EOS of Fraga, Pisarski and Schaffner-Bielich and our EOSs begin to 
show appreciable difference: the pressure density given in our EOSs is much lower than 
that given in the EOS of Fraga, Pisarski and Schaffner-Biehch. Here we note that the 
form of our EOS depends strongly on the model quark propagator at zero /i one chooses. 
The parameters of the model quark propagator (9) are determined by numerical fitting [25] 
and there exist some arbitrariness in this process. In further researches we will apply our 
EOS to the study of neutron star and wc hope that by comparing theoretical results with 
observational results of neutron star we can have further constraints to the parameters of 
our model quark propagator. 

To summarize, in the study of strongly interacting matter at finite temperature and/or 
finite density, the knowledge of the partition function of QCD determines all the thermody- 
namic properties of the system. In this paper, we try to give a direct method for calculating 
the partition function, and hence the EOS of QCD at finite chemical potential and zero 
temperature. In this method the quark-number density p(/i) is expressed in terms of the 
dressed quark propagator at finite chemical potential G[ji]{p), and the pressure density V{ii) 
is given by the integration of p{ii) plus an additive constant, whose physical meaning is the 
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FIG. 4: The pressure, relative to the free quark gas pressure Vfree = -^c^//^^/(12vr^), in our EOS 
(18) and the EOS (22) of Fraga, Pisarski and Schaffner-Biehch 

pressure density at zero chemical potential 'P(/i)|^=o. Because of the difficulty of calculating 
G[ij]{p) and V{fi)\^=o from first principles of QCD, one has to resort to nonperturbative 
QCD models when applying our method. In this paper we adopt one nonperturbative QCD 
model, the rainbow-ladder approximation of the Dyson-Schwinger approach. In order to 
obtain an EOS with explicit analytical form, we choose the meromorphic model quark prop- 
agator proposed in Ref. [25]. Using the general result proved in the framework of the 
rainbow-ladder approximation of the DS approach in Refs. [17,18], G[fi]{p) is obtained from 
this model quark propagator. From this the quark-number density p(/i) is calculated, which 
is found to differ significantly from the Fermi distribution of free quark theory. Physically 
this is a consequence of dynamical chiral symmetry breaking and confinement in the low 
energy region. It is found that when fi is below a critical value /io (yUo = 351 MeV, 377 MeV 
and 341 MeV for the 2CC, IRICC and 3R parametrization, respectively), the quark-number 
density vanishes identically. This feature agrees with the general conclusion in Ref. [27]. 
The value /iq obtained here is almost of the same order of magnitude as the estimate made in 
[27] (307 MeV). The constant P(/i)|^=o is also self-consistently calculated using the rainbow- 
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ladder approximation of the DS approach. Prom these the full analytic expression of the 
EOS of QCD at finite and zero T is obtained (apart from the constant term P(//)|^=o 
which can in principle be calculated from the C JT effective action) . A comparison between 
our EOS and one prominent example of the EOS of QCD, the cold, perturbative EOS of 
QCD proposed by Fraga, Pisarski and Schaffner-Bielich in [30] is made. 
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